ASYMPTOTIC BEHAVIOR OF SOLUTIONS TO THE LIQUID 
CRYSTAL SYSTEMS IN H m (R 3 ) 

MIMI DAI AND MARIA SCHONBEK 

Abstract. In this paper wo study the large time behavior of regular solutions 
to a nematic liquid crystals system in Sobolev spaces i? m (R 3 ) for m > O.We 
obtain optimal decay rates in H m (R 3 ) spaces, in the sense that the rates 
coincide with the rates of the underlying linear counterpart. The fluid under 
consideration has constant density and small initial data. 



1. Introduction 

In this paper we consider the asymptotic behavior of solutions to the simplified 
model of nematic liquid crystals (LCD) with constant density in Sobolev spaces 
H m (R 3 ) for m > 0: 



u t + u ■ Vu + V7r = isAu - V • (Vd (8 Vd), 
(1.1) d t + u-Vd = Ad-f(d), 

V • u = 0. 

The equations are considered in R 3 x (0,T). Here, n : R 3 x [0, T] — > R is the fluid 
pressure, u : R 3 x [0, T] -> R 3 the fluid velocity, and d : R 3 x [0, T] -> R 3 the director 
field representing the alignment of the molecules. The constant v > stands for 
the viscosity coefficient. Without loss of generality, by scaling, we set v = 1. The 
forcing term Vd <g> Vd in the equation of conservation of momentum denotes the 
3x3 matrix whose ij-th entry is given by "V;d ■ Vjd", 1 < i, j • < 3. This force 
Vd (g> Vd is the stress tensor of the energy about the director field d, where the 
energy is given by: 



|Vd| 2 da;+ / F(d)dx 



where 



F(d) = ^(|d| 2 - l) 2 , f(d) = VF(d) = ^(|d| 2 - l)d, 



for a constant rj in this paper. The F(d) is the penalty term of the Ginzburg-Landau 
approximation of the original free energy for the director field with unit length. 
We consider the following initial conditions: 

(1.2) u(x,0) = u (x), V-ti = 0, 

(1.3) d(x,0) = d (x), |d (x)| = l, 
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and 



(1.4) 



u e H m (R 3 ) 7 d - w e H m+1 (R 3 



) 



for any integer m > 1 with a fixed vector wq ^ S 2 , that is, |u>o| = 1. 

The flow of nematic liquid crystals can be treated as slow moving particles where 
the fluid velocity and the alignment of the particles influence each other. The hy- 
drodynamic theory of liquid crystals was established by Ericksen [8j and Leslie 
P13 [16] in the 1960's. As F.M. Leslie points out in his 1968 paper: "liquid crys- 
tals are states of matter which are capable of flow, and in which the molecular 
arrangements give rise to a preferred direction" . There is a vast literature on the 
hydrodynamic of liquid crystal systems. For background we list a few, with no 
intention to be complete: H21 [13l [HI HH HT1 [U EJ [3j [H EU [20l [H]. In particu- 
lar, the asymptotic behavior of solutions to the flow of nematic liquid crystals was 
studied for bounded domains in |18l HM] . It was shown in that, with suitable 
initial conditions, the velocity converges to zero and the direction field converges 
to the steady solution to the following equation 



In [24] . Lemma 2.1 the Lojasicwicz- Simon inequality is used to derive the conver- 
gence when 51 is a bounded domain. 

In our previous work [5], we established a preliminary decay rate for the solutions 
in K 3 to (|1.1[) , subject to the additional condition on the director field which insures 
that the initial director field tends to a constant unit vector wq, as the space variable 
tends to infinity: 



This behavior at infinity of the initial director field allows to obtain the stability 
without needing the Liapunov reduction and Lojasiewicz- Simon inequality, since 
Wo is a non-degenerate steady solution to (|1.5[) . 

The paper is organized as follows, in section 2 we recall some previous results, 
give some preliminary estimates and state the main Theorem. In section 3 we ob- 
tain the decay rate for the velocity in L 2 (M. 3 ). Section 4 deals with the decay rate 
for gradient of the director field in L 2 (M. 3 ). The last section gives decay rates in 
H m for the velocity and the director vector. All the obtained rates are optimal. 
The rates obtained improve the rates obtained in [5]. We work with regular solu- 
tions with small data. The main tools used are the Fourier splitting method and 
appropriate energy estimates. 

Remark 1.1. We note that for the L 2 the decay rates can be obtained also for 
weak solutions. In this case one shows the decay for approximations in the form we 
do in the paper and passing to the limit it will follow for the solutions. We expect 
that for higher derivatives one can obtain the decay for approximations, and show 
that eventually the solutions become small and the use the results in this paper. 



(1.5) 




(1.6) 



lim do(x) = Wq. 

| x I — >■ oc 
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2. Preliminaries 

We start from the basic energy estimates and Ladyzhcnskaya higher order en- 
ergy estimates [HI [5] . We establish that the velocity u converges to with decay 
rate (1 + t) -1 / 4 in L 2 (R 3 ) and that the convergence of the direction field d to the 
constant steady solution wq with the decay rate (1 + t)~' s ^ 1 ~^ in L P (R 3 ) for any 
p > 1. The main ingredients are the Fourier splitting method and the Gagliardo- 
Nirenberg interpolation techniques. 



We recall that for small data there exists a smooth solution given by 

Theorem 2.1. [5] Let u and do satisfy U.2\) - [L4\ l. Assume that uo G ff 1 (R 3 ) 
and do — Wo G H 2 (R 3 ) n L 1 (K 3 ) for a unit vector w . There is a positive small 
number eo such that if 

( 2 -7) IM|ffi( R 3) + IMo ~ wo\\ 2 H 2 {R3) < eo, 

then the system il.l]) has a classical solution (u,Tr,d) in the time period (0,T), for 
all T > 0. That is, for some a G (0, 1) 

u G C 1+a/2 ' 2+a ((0,T) x R 3 ) 

(2.8) Vtt G C a/2 ' a {{0,T) x R 3 ) 

deC 1+a/2 - 2+a ((0,T) xR 3 ). 

And the solution (u, tt, d) satisfies the following basic energy estimate and higher 
order energy estimate 

(2.9) / \u\ 2 + \Vd\ 2 + 2F(d)dx + 2 f f |Vu| 2 + \Ad - f{d)\ 2 dxdt 

< INH^RS) + ||Vdo||i2( R 3) 

(2.10) / \Vu\ 2 + \Ad\ 2 dx + [ [ \Au\ 2 + \\7Ad\ 2 dxdt 
Jr 3 Jo jr3 

< C(I|W0|Ih1(R3) + |M0 - ^0|Hf2 (H 3)). 

with the constant C depending only on initial data and on r\. 

In [5] we established the following decay result for the regular solutions obtained 
in Theorem O 

Theorem 2.2. Let (u,n,d) be a regular solution to the system M.l\) with initial 
data \1.2)) and \1.3\) and boundary condition {1.4]) for m = 1. Assume additionally 
uo G L X (R 3 ) and do — wo G L P (R. 3 ), for any 1 < p < oo and a unit vector wo- There 
exists a small number eo > such that if 

(2-11) \\ u o\\m(K3) + \\do ~ wo II 2(1*3) < e , 

then 

(2.12) K-^)-™o||lp(R3) < C p ||d -Wo||LP(R3)(l + t)"^ 1 -5), 

(2.13) \\V(d(-,t)-wo)\\h (R3) <C(l + t)--t, 
(2-14) ll«M)lli»™<C(l +*)"*, 
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where the various constants C depend only on initial data, C p depends on the data 
and p. 

We first need to establish energy estimates for the velocity u in H m (M. 3 ) and the 
director field d in H m+1 (M. 3 ), for all m > 1, starting with more regular initial data. 
We will use the following notation 

(2-15) *l(t) = \\D k u\\ 2 L2{R3) + \\D k+1 d\\ 2 L2{m . 

m 

(2-16) *™(*) 

k=0 

Theorem 2.3. Let (u,7r,d) be a regular solution to the system M.l)) obtained in 
Theorem \2.1[ Assume {l-4ty holds. Then the solution (u, n, d) satisfies, for all 
m > 1 

(2-17) / V 2 m dx+ [ [ ^ 2 m+1 dxdt 

< C m (\\u \\ 2 Hm(R3) + 1 1 d - w \\ 2 Hm +i (R3) ). 

The constant C m depends only on initial data, rj and Ck for k = 0, . . . , m — 1, and 
Co depends only on the data. 

Our main result reads as follows 

Theorem 2.4. Assume the initial data (uo,do) satisfies n.2\ )- [T7^) , Uo 6 i 1 (M 3 ) 
and d — wq g L p (R 3 ) for p > 1. Let (u,TV,d) be the regular solution obtained in 
Theorem \2.1\ Then, for all m > 



(2.18) \H;t)\\ 2 Hm{R3) <c m (i + t)-( m +i\ 

(2.19) \\d(;t)-w \\ 2 Hmm <C m (l+t)-^+i\ 

(2.20) \\D m u(;t)\\ L ~ {R3) <C m (l + t)-^, 
and 

(2.21) \\D m (d(;t)-w )\\ L oo m <C TO (l+t)-^. 



Here the various constants C m depend only on initial data and m. 

Remark 2.5. In [5] where we established Theorem l2.2[ we remarked that inequality 
(|2.12|) docst not include the limit case p = oo which would produces the decay of 
d — uiQ in i°°(K. 3 ) with the rate (1 +i)~ 3 / 2 . The reason was that the small constant 
eo in Theorem 12.21 would be forced to be zero if p = oo. In this paper, we are able 
to obtain the optimal decay rate (1 +t)~ 3 / 2 for d — wo in L°°(R 3 ), which is included 
in (12311 . 

In [21] , the authors studied the large time behavior of solutions to Navier-Stokes 
equation in H m (M. n ) for all n < 5. The decay estimate (|2.18[) for the velocity 
obtained in Theorem l2.4l coincides with the result in [21] for Navier-Stokes equation. 



ASYMPTOTIC BEHAVIOR OF LCD IN H m (R 3 ) 



5 



3. Higher order energy estimates 

In this section we show that solutions with initial data in H m (M. 3 ) remain in 
H m (R 3 ) for all time. We prove the energy inequality ([2~T7|) in Theorem [231 The 
main tool used is a modified Ladyzhenskaya energy method (see [6]). 

In the sequel we need to use a Gagliardo-Nirenberg interpolation inequality. For 
completeness we recall the inequality here 

Proposition 3.1. Q~U] Let w G W m,p (K n )nL 9 (K n ), for 1 < p < oo and 1 < q < oo. 

Then 

(3-22) \\D k w\\ Lr{Rn) < C^D^l^Ml-^ 

for any integer k S [0, m — 1], where 

1 k , 1 m s , ,1 

3.23 - = - + a( ) + (i_ a )_ 

r n p n q 

with a € 1], either if p = 1 or p > 1 and m — k—- J\fU {0}, while a <E 1), 
ifp > 1 andm- k-% e Afu {0}. 

We now prove Theorem 12.31 

Proof of Theorem 12.31 The proof proceeds by induction. We recall that 
inequality (|2.17p for m = is easily obtained by energy estimates and has been 
established in [5]. Suppose (|2.17|) holds for m = 1, 2, ... ,j — 1, that is 

(3.24) / ® 2 - 1 dx + [ [ ^dxdt 

Jw Jo Jr 3 

< Cj_l(||Uo||ff m ( ffi 3) + \\d - W0||ffm + 1( K 3)). 

We need show (f2~T71) holds for ?n = j. Let < fc < j. Recall that 



(3.25) & k (t) = ||D fc U ||i 2(H3) + \\D k+1 d\\ 2 



The idea is to establish an ordinary differential inequality for $|(f) usm g the 
equations in system and Gagliardo-Nirenberg interpolation inequalities, and 

then sum all the terms <&§(i), $f (i), . . . , $\{t). 



Remark 3.2. We note that since u, d are in H m and are regular, there are no 
boundaries that we have worry about in the subsequent integration by parts. 

S2 



Taking derivative of (t) with respect to time yields 



dx 



(3.26) -^-<f> 2 k (t)= [ D k u-D k u t dx+[ D k+1 d- D k+1 d t dx 

2 dt J R 3 J R 3 

= -[ D k+1 u- D^utdx- ( D k+2 d-D k d t 

JR3 JR3 

Taking the (k — l)-th spatial derivative on the first equation in gives 
D k ~ l ut = D k ~ 1 Au - D k ~ 1 (u ■ Vu) - D fe-1 V7r - D* _X (V • (Vd® Vd)). 
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Thus, 
(3.27) 



D k+1 uD k ~ l u t 



= - [ \D k+1 u\ 2 dx+ [ D^uD^iu-V^dx 
Jr 3 Jr 3 

+ f D k+1 uD k - 1 {V -{Vd®Vd))dx 
Jr 3 

<— I \D k+1 u\ 2 dx + A I \D k (u®u)\ 2 dx + 4 I \D k {Vd®Vd)\ 2 dx. 
2 Jr 3 Jr 3 Jr 3 

Taking the fc-th spatial derivative on the second equation in gives 
D k d t = D k Ad -D k [u- Vd) - D k {f{d)). 

Thus, 
(3.28) 

D k+2 dD k d t 

{ \D k+2 d\ 2 dx+ [ D k+2 dD k {u-Vd)dx + [ D k+2 dD k f{d)dx 
Jr 3 Jr 3 Jr 3 

-f \D k+2 d\ 2 dx + A [ \D k {u-Vd)\ 2 dx + A\ f \D k ((\d{ 2 -l)d)\- 
2 Jr 3 Jr 3 V J R 3 

Inserting (f3T27|) and ff3T28|) into (f3~26]) yields, 

(3.29) x $2 W+ / \D k+1 u\ 2 dx+ [ \D k+2 d\ 2 dx 

at J K 3 J R3 

\D k (u®u)\ 2 dx+ ( \D k (\7d®S7d)\ 2 dx 



< 



dx. 



< 4 



+ 4 



f \D k {u-Vd)\ 2 dx+\ f \D k ((\d\ 2 -l)d)\ 2 dx 
Jr 3 V Jr 3 



= 4[7 + II + III + IV}. 

The four terms on the right hand side of (|3.29[) are estimated as follows. 
For /, we have 



(3.30) 
with 

(3.31) 



J<2H=Lo.||i?*ii||i a +2||Vu||i.||I>*- i «||i a +7 1 



fc-l„.l|2 



0, if fe = 2,3 

Ei ( =2 2)721 IP r «lli-IP fc_r « 



k ~ r " A %, if fe > 4. 



Where [a] denotes the largest integer that is less than or equal a. 

Recall u is regular by (|2.8p . ||m||z,°° and || V«||l°° are bounded by some fixed constant 

C. By the induction hypothesis (pT24"]) . ||L»' £ " 1 it||2 2 is bounded by Thus, 

(3.32) ||«||i-||i? fc «||i a + HV-ulliooll^-^Hi, < C[*|(t)]. 

By the remark above C denotes an absolute constant. Applying Gagliardo-Nircnberg 
interpolation inequality, 

\\D r u\\ L ~ < C\\D r+2 u\\h\H^ a 
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with a = j-Q < 1. Thus, if k > 4, by the last inequality and p. 311) . 



2r+4 

p-2)/2] 

(3.33) h<C J2 \\D r+ Ml^D k - r u\\l 2 \\uC- 1 

r=2 

Note that when k > 4. r + 2 < fc — 1 for all r < [(k — 2)/2]. Hence, by induction 
hypothesis (|3.24p . it follows from (|3.33[) 

(3.34) h < C k ^l(t). 
Combining (f3"3Djl . ([331) . (£52)1 and (jXM)) gives 

(3.35) I < C k * 2 k (t) 
For /J, we have 

(3.36) II < 2||Vd|||oc \\D k+1 d\\ 2 L2 + 2\\D 2 d\\ 2 L ^ \\D k d\\\ 2 + 77 1; 
with 



(3.37) 77 x = 



0, if = 2, 3 

El=2 2)/2] IP r+1 rf|l!oop fe - r+1 d|li 2! if ft > 4. 



(3.39) //i<C P r+3d ll!"P fc_r+ld lll=ll Vd llj (1 



Since d is regular in the sense of (|2.8j) . ||Vd||L°° and |7? 2 d||£°° are bounded. By 
induction hypothesis (|3.24|) . |7> fc <i|||2 is bounded by Thus, 

(3.38) ||Vd||i« \\D k+1 df L , + \\D 2 d\\U \\D k d\\l 2 < C*l(t). 

Applying Gagliardo-Nirenberg interpolation inequality, 

\\D r+1 d\\ L ~ < C\\D r+3 d\\ a L 4Vdf£- a) 
with a = |i±| < 1. Thus, if k > 4, by the last inequality and (|3. 371) . 

[(fc-2)/2] 

II n'-+3^l|2a || n fe-r+l^||2 ||r7JI 
r=2 

Note that when fe > 4, r + 3 < ft for all r < [(ft — 2)/2]. Hence, by induction 
hypothesis (|3~2"4"]) . it follows from (|3~39l 

(3.40) 77 x < C fc ^(t) 
Combining (f3~36l> . ([33l|), (j3~38| and (j3~40| gives 

(3.41) 7/ < CfctfjUt). 
For 7/7 , we have 

(3.42) 777 < ||u||^||7) fc+1 d||l 2 + || Vu||£« ||7? fc d|l^ + P 2 "lli~P^II^ 

+ llVdlllooll^Hi, + ||D a d||i 00 ||2?*- 1 u||i a +777 1; 

with 

0, if ft = 2, 3, 4, 

J2r=l\\D r d\\U\D k -^u\\l 2 , if fc>5. 

Recall u and d are regular by (|2.8|) . ||7) r 'u||i,oo and ||7? r d||Lo° are bounded for r = 
0,1,2. By induction hypothesis ([3T2Ij) . ||7> fe - 1 u||| 2 , || J D fe " 1 d||| 2 and ||7> fe d||| 2 are 
bounded by Thus, by fl3H]), 

(3.44) 777 < C^l(t) +777i. 



(3.43) 777 x 
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Applying Gagliardo-Nirenberg interpolation inequality, 

(3.45) \\D r d\\ L - < C\\D r+2 d\\l 2 \\d - woW 1 "^ 
with a = f£±§ < 1. Thus, by (15~33l 

(3.46) m^C^V^II^P^^II^IIrf-^ollii 1 -^. 

r=3 

Note that r + 2<k and k - r + 1 < k - 2 for all r G [3, k - 2], when fc > 5. Hence, 
by induction hypothesis (|3.24j) . it follows from ()3.46j) 

(3.47) II h <C k *l. 
Combining flgTUD and (|3~4"T)) gives 

(3.48) III <C$ 2 k (t) + C k . 
Note that |d| < 1. For IV, we have 

(3.49) IV < C\\D k d\\\ 2 + llVdllloclp*-^!!^ + \\D 2 d\\l x \\D k - 2 d\\l 2 

k-3 

+j2\\ Drd \\i°°\\ Dk - rd \\h- 

r=3 

By (f3T45|) and induction hypothesis (|3~24| . it follows from ([349]) 

(3.50) /V < C fe ^(i). 
Combining (j3~2^)) . (j3"35)) . (j3~4Tj) . (|3~4^|) and (j3~5(J)) yields 

|$ 2 fe (t) + \\D k+1 u\\l 2 + \\D k+2 d\\ 2 L2 < C k *l(t). 
Summing the last inequalities for k = 0, . . . , j yields 

(3.51) ±^ j (t) + ^ +1 (t)<C j ^ j (t). 

Integrating (|3.5ip yields by the inductive Hypothesis 
(3.52) 

/ ^ *K°) + C 'j / ^ C J (II«oIIh*(R») + No - W0||Hi+i(R3)) ■ 

Jo Jo 
This concludes the proof of the Theorem. 

□ 



4. Improvement of the velocity decay rate 

In [5], for solutions with uo G L 2 (R 3 ) flL 1 (R 3 ) , appropriate assumptions on the 
data do and sufficiently small data, we showed that the velocity satisfies ||it|| 2 2 < 
C(l + In this section, with the same assumptions we show that the velocity 

has the optimal decay rate ||u|| 2 2 < C(l + t)^^ . Namely, we show 

Lemma 4.1. Assume u a G H^R 3 ) n ^(R 3 ) and d - w G H 2 (R 3 ) n L^R 3 ). Let 
u be a smooth solution obtained in Theorem \2.1\ with data satisfying \2.11\) . Then, 
u satisfies 

(4.53) Nl! 2(R 3)<C(l + i)-i. 
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Proof: As in [5], the Fourier splitting method, [22l[23], is applied to establish the 
decay of solution. We decompose the frequency domain R 3 into two time dependent 
subdomains S(t) and its complement S c (t). Define 

(4.54) S(t) = U G K 3 : 1^1 < r(t) = { — ^ 2 



i + r 

for a constant k that will be specified bellow. One of the key estimates to establish 
decay of velocity 

(4.55) Mi»( B 3)<C(l + *)"* 

in [5] is that the Fourier transform of the solution u satisfies 

\u(tt)\ < Cltr 1 for£eS(t). 

We show first that with the decay estimate (|4.55j) . the Fourier transform of the 
velocity u satisfies 

(4.56) \u(£,t)\<C, for^eS(i) 

for an absolute constant C . As a consequence, we proceed with a similar analysis 

as in [5] and obtain the optimal decay (|4.53p for the velocity. 

Taking the Fourier transform of Navier-Stokes equation in system (jl.ip yields 

(4.57) u t + \^\ 2 u = G(^t) 
where 

G(£, t) = -T(u • Vu) - ^(Vtt) - J-(V • (Vd ® Vd)), 
and J 7 indicates the Fourier transform. Multiplying ()4.57|) by the integrating factor 
e l£l * yields 

|[ e l^] = el^G(£,i). 

Integrating in time gives 

/•t 

(4.58) u(f,i) = e - | « |2t w + / e-^ 2{t - s) G{^s)ds. 

Jo 

We analyze each term in G(£, <) separately. We have 

(4.59) |.F(u- Vu)| = |J"(V- («®«))| < V / |«V||^-|dx < C(l +t)~ 1/2 |C| 
due to (|435l> . 

In [5], we obtained decay for the gradient of director field d as 

||Vd||i 2(R3) <G(l+t)" 3 / 4 . 

Thus, we have 

(4.60) | J"(V • (Vd ® Vd))| < C(l + i)~ 3/4 |£|. 
Taking divergence of the velocity equation in system (jTTTJ) yields 

Taking the Fourier transform of the last equation gives 
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Combining ()4.59|) . (|4.60j) and the last equation yields 

JT(7r) < C(l + t)" 1 / 2 , 

and thus 

(4.61) T(Wtt) < C(l +t)- 1 / 2 \£\. 
Combining (|4~59]) . (|4T60j) and (|4~6Tj) yields 

(4.62) < C(l + t)- 1/2 |^|, for£eS(t). 
From (|4.58p and (|4.62[) we have 

(4.63) |u(£,t)| <e-'« |2t K|+ f e-^ 2it - s \l + S r^ 2 \^\d S . 



Since uq G L , we have |uo| < C for all £. Performing integration on the right hand 
side of (|4.63p gives 

m,t)\ < Ce-I«l 2 ' + C(l + i) 1/2 |CI < C 
since £ G and |£| < (7(1 + 1)" 1 / 2 . This completes the proof of (OB) , 

Multiplying the velocity equation in (jl.ip by u and integrating by parts yields 
1 d 



2dt 



/ |w| 2 d.T+ / \Vu\ 2 dx = / \7u(Vd®Vd)dx 

JR 3 Jr 3 JR 3 

1 / |Vu| 2 dx + C / \\7d®Vd\ 2 dx. 

2 7r 3 Jr 3 



< 

Thus, 

d , 2 



/ \\7u\ 2 dx<c[ \Vd®Vd\ 2 dx. 

JR 3 Jr 3 

The right hand side of (|4.64[) is estimated by 



(4.64) -/ |„ 



/ \Vd(g> Vd\ 2 dx = / (\7d(g)Vd){Wd<E)Vd)dx 
Jr 3 Jr 3 

= -3/ (d-w ) AdVd®Vd 
Jr 3 



< 

It follows that 



\ I \\7d®\7d\ 2 dx + C [ \d- w \ 2 \Ad\ 2 dx. 
2 7r 3 7r 3 



/ \Vd®Vd\ 2 dx<C [ \d- w \ 2 \Ad\ 2 dx 
Jr 3 Jr 3 

<c(J \d~w \ p dxy (J \Ad\^dxj 

= C\\d-w \\ 2 LPm ( [ \Ad\ 2+ ^dx 
\Jr 3 



< C\\d-w \\ LP{R3) 
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for p > 2. The last step followed from the energy estimate (|2.10j) and recalling that 
||Ad||ioo(R3 X [ 0l T]) is bounded since d G H m (R 3 ) is regular by (|2.8j) in Theorem |2~T1 
Thus, ([2T2]) in Theorem \T2\ yields 

/ |Vd® X7d\ 2 dx < C(l + t)- 3{1 -p\ 
Jr 3 

for any p > 2. Therefore, it follows from (|4.64l) 



(4.65) 



^- ( \u\ 2 dx+ [ \Vu\ 2 dx <C(l + t) 

dt J K 3 J R 3 



-3(l-i) 



Applying Plancherel's theorem to (|4.65j) gives 



/ n 2 di+ ( iei 2 N 2 rfe<c(i+*)" 3(1_ ^- 

JR 3 JR 3 



(ft 

We reorganize the last inequality as 

(4.66) ^/ \u\ 2 dt<-[ \e\u\ 2 d^- [ \e\u\ 2 dt + C(l + t)- 3(1 -^ 
dt Jr 3 Js(t) c JS(t) 

k 



< -- 



1 + i 



\u\ 2 d£- / |^| 2 |M| 2 dC + C(l+0" 3(1 ~^ 



<-Az/ N 2 ^ + t^/ H 2 ^ + c(i + i)- 3(1_i) - 

S(t) 



1 + i 7 R 3 1 + f 

Due to (|4.56[) we have |w| < C, hence 

,r(t) 



/ \u\ 2 d£<C r 2 dr < C(l +t)~ 3/2 . 

.Asm 



By (|4.66[) choosing any p > 6 



dt J R 3 1 + t 

< G{l + t)-i. 
Multiplycation by the integrating factor (1 + t) k yields 
d 



(1+t)* \u\ 2 dt 



< C{l + tf 



dt 

Integration in time gives 

(l+t) k [ \u\ 2 d£< I \u(£,0)\ 2 dZ + C[(l + t) k -i -1}. 



Thus, 

f \u\ 2 dt<(l + t)- k f \u^0)\ 2 d^ + C[(l + t)-l -{l + t)- k }. 
jr 3 Jr 3 

Since u G L 2 , u(0) G L 2 by Plancherel's theorem. We choose k > 3/2. Hence 

\u\ 2 dx < C{l + t)-i. 
The proof of the Lemma is complete. 

□ 

The rest of this section deals with an auxiliary estimate for the velocity u in L°°. 
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Lemma 4.2. Let u be the solution of system 1 1.1]) with initial data satisfying the 
conditions in Theorem \2.1\ Then we have 

(4.67) \\u\\ Laom <C m {l + t)-^ l -^\ 

for m > 1 . 

Proof: By Gagliardo-Nirenberg inequality p.22[) . we have 

IMU- < \\D m u\\i4u\\ l -s 

with a = Combining the estimates (|2.17j) . (|4.53|) and the above inequality 

yields 

HU- <c TO (i + f)-* (1 -^ 

which establishes the conclusion of the lemma. 

□ 

Remark 4.3. The rate in Lemma |4~21 will be improved to (1 + t)~ 3 > 2 in Section^ 



5. Optimal decay rate of Wd 
The goal of this section is to establish that Vd decays in L 2 (R 3 ) at the rate 

(i + t)-i. 

We first obtain an auxiliary decay rate for d — wq in i°°(R 3 ), 

Lemma 5.1. Let d be the solution obtained in Theorem \ 2.1{ Assume the initial 
data satisfies the conditions in Theorem \2.2i Then 

(5.68) \\d(-,t)-wo\\L~p*)<C(l + t)-i-%*, p>l 
where the constant C depends only on the initial data. 

Proof: By the Gagliardo-Nirenberg interpolation inequality, 

Hd-ttfoiu- < ei£ 3 d||M|d-™ollL; a 

with a = 2^ ■ Due to the energy estimate (|2.10|) and the decay estimate in Theorem 
12.21 we have that 

\\d - w ||l~ < C(l + t)-f C 1 "!)^-^), p > 1, 
which proves (|5.68|) . 

□ 

Corollary 5.2. Let d be the solution obtained in Theorem l2.il Assume the initial 
data satisfies the conditions in Theorem 12.21 Then 

(5.69) (d + w ) ■ d>0, 
for sufficiently large time t. 

Proof: Taking p = 7 in (|5.68[) yields 

\\d-wo\\ L oo KCil+t)- 1 . 
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Hence, we have 

(d — wo) ■ d = \d\ 2 + d ■ wo 

= \d\ 2 + (d~w )-w + l 

> \d\ 2 + 1 - \d - w \\w \ 

> \d\ 2 + l-C(l+t)- 1 

> 

for sufficiently large time t. This proves (|5.69l) . 

□ 

Lemma 5.3. Let d be the solution obtained in Theorem \ 2.1[ Assume the initial 
data satisfies the conditions in Theorem \2.2i Then 

(5.70) \\Vd(;t)\\ L2{m <C(l+t)-i. 

The constant C depends only on the initial data. 

Proof: Multiply the director field equation in by Ac? and integrate over 

M 3 , then 

( 5J1 ) \i I Nd\ 2 dx+ [ \Ad\ 2 dx = [ {u-Vd)Addx + [ f(d)Addx. 

2 at J R 3 J R 3 J R 3 J R 3 

We need the following auxiliary estimates 

/ (it • Vd)Addx < - [ \Ad\ 2 dx + C\\u\\ 2 L ^ [ \Vd\ 2 dx 

JM? 4 J M 3 J R 3 

and since \wq\ = 1 



/ f(d)Addx = -\ [ V[(d + w )(d - w )d]\7ddx 

Js.3 If JW 



< C\\d-w \ 



L- 



[ \Vd\ 2 dx-\( (d + w )-d\\7d\ 2 dx 
7r3 r\ J R 3 



< C\\d-w \\L°° [ \Vd\ 2 dx, 



where we used the facts that \d\ < 1 and (d+wo) -d > for large time, see Corollary 
5.21 Combining the last two inequalities with (|5.71j) yields 



— / \Vd\ 2 dx+ [ \Ad\ 2 dx 

dt Jg3 J R 3 

< C(l|«|lioo + \\d-w \\ L °o) I \Vd\ 2 dx 



KCiil + t)-^ 1 -^ + (l + tyi-^) I \Vd\ 2 dx 



due to Lemma 14.21 and Lemma |5. II Taking m = 5,p = 7 yields that 

f \Ad\ 2 dx< Coil + t)- 1 [ \Vd\ 2 dx. 

JR3 JR3 



(5.72) j t I 



Apply now the Fourier splitting method to obtain the decay estimate. We choose 
an appropriate constant k for the time dependent sphere S(t) as in (|4.54|) . such 
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that k — Co > 5/2. Here the constant Co is the one on the right hand side of (|5.72[) . 
Proceeding by the Fourier splitting method gives 



(5.73) 



/ \Ad\ 2 dx> f \e\T(vd)\ 2 dt 

Jr 3 Jr 3 \s 

> -A- / \F(Vd)\*d£ 



A- 



1 + i 

Combining (f5~72"l) and (|5~73| yields, 



/ |J-(Vd)| 2 ^--^- / |J-(Vd)| 2 ^. 



d /" ,|2 . , fc _ 



(5.74) — / \\7d\ 2 dx+- -/ \Vd\ 2 dx 

Ut J R 3 1 + t J R 3 

Multiplying (|5.74p by the factor (1 + t) k ~ Co and integrating in time yields 

/ \Vd\ 2 dx < (l + t)" (fe " Co) / iVdo^dx + il+ty 1 [ \d-w Q \ 2 dx 
Jr 3 Jr 3 Jr 3 

<C(l+t)-i, 

due to the decay estimate in Theorem 12.21 and the choice k — Co > 5/2. This 
completes the proof of the Lemma. 

□ 



6. Decay of solutions in higher order Sobolev spaces 

In this section we obtain the decay estimates for u in H m (M. 3 ) and d in H rn (M. 3 ) 
with m > 1. The method involves induction and a Fourier splitting argument. 

We recall the following Gagliardo-Nirenberg inequalities 

(6.75) W&uWloo < CWD^uW^hW 1 -^ 

(6.76) H^dlU- <C||D ro+1 d||^||d-«;o||^ 
with a. t = for i > 1 and 1 + 3/2 < m - 1/2. . 

We first establish the following auxiliary estimates. 
Lemma 6.1. Let d be the solution obtained in Theorem \2.1{ Then d satisfies 

(6.77) ||Vd||i°o < C(l + i)-t(i-A) 

(6.78) \\D 2 d\\ L °o < C(l + t)-i^~A) 
for k > 1. 
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Proof: Taking w = Vd, fc = 0, r = oo, and p = 9 = 2 in (|3.22[) yields, for fc > 2 

||vd|| L » <cp fe rf|IMIvrf|ll7 a 
<c(i + t)-3 (1 -* ) 

due to the higher order energy estimate (|2.10[) in Theorem 12.11 and the decay es- 
timate (|5.70[) in Lemma T5.3I The constant C depends only on initial data. Thus 
(|6.77p is proved. The proof of (|6.78p is similar and as such is omitted. 

□ 

Remark 6.2. The decay of Vd and Ad in i°°(IR 3 ) will be improved to the optimal 
rates. 

Wc now establish a higher order energy estimate for the solution. The ideas are 
based on work in |21j . 

Lemma 6.3. Let (u,d) be the solution obtained in Theorem ] 2. 11 Assume the initial 
data satisfies the conditions in Theorem \2.2[ Then (u,d) satisfies, for m > 1 



(6.79) f \D m u\ 2 + \D m d\ 2 dx+ [ \D m+1 u\ 2 + \D m+1 d\ 2 dx 

dt J R 3 J R 3 

< C m (Nli- + \\d- «*>||i- + \\D 2 d\\ 2 L ^) f \D m u\ 2 + \D m d\ 2 dx 

JTR3 



CmWd-WoWlcoWD^dWlz+Rr, 



where 
(6.80) 

'0, m = l,2; 



Rn 



ESi MUiD^uWir* lld-woiiUD^-^df-^ 
+ TZ ~i \\d - M\h \\D m -*u\\Jr + e2t ¥ - Ml* h^-miits 

^ = R m l + Rm2 + R,n3 + #ro4, for 771 > 3 

with a,i = ra+i ■ ^ e constants C in the inequality depend only on the initial data. 

Proof: Taking the m-th derivative on the first equation in multiplying it 

by D m u and integrating over R 3 yields 



1 d 
2dt 



D m u\ 2 dx + I \D m+1 u\ 2 dx 



< 



- [ D m (u ■ Vu)D m udx - [ D m {V ■ (Vd ® \7d))D m udx 

JR3 JR 3 

/ D m - 1 (u-\7u)D m+1 udx+ [ D m {Vd®\/d)D m+1 udx 

JR3 JR3 

- [ \D rn+1 u\ 2 dx+ [ \D rn - 1 (u-\7u)\ 2 dx+ [ \D m {Vd®Vd)\ 2 dx. 

4 JR3 Jr3 7 H 3 
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Therefore, 

(6.81) ^ / \D m u\ 2 dx + ^ [ \D m+1 u\ 2 dx 
dt JR 3 2 J R a 

< [ \D m - 1 (u-\7u)\ 2 dx+ [ \D m {Vd®Vd)\ 2 dx 
Js 3 Jr 3 

< ll«||i-||D m «||i a + llVdlli-p^dlH, + \\D 2 d\\ 2 ^\\D m d\\ 2 L2 

+ E 11^111.11^-^1113+ \\D l+1 d\\U\D m -^d\\ 2 L2 . 

l<i<m/2 2<i<(ra-2)/2 

Using (|6.75l) gives that 

(6.82) J2 ll^lll^p—^Hl, 

l<i<m/2 

l<i<m/2 

<\\\D m+1 u\\ 2 L2 + C Y, MUD™-^ 1 -^. 

l<i<m/2 

Using (|6.76p gives that 



(6.83) E ll^dlli-ll^^dlli- 

2<i<(m-2)/2 

<c E \\D m+1 d\\ 2 £ +i \\d-w4 2 £- a ^\\D m - i +H\\ 2 L2 

2<i<(m-2)/2 

< \\\D m+1 d\\ 2 L . +C £ lld-^H^P"- 4 + 1 d\\$ 1 -»* ) . 

2<i<(m~2)/2 

Taking the m-th derivative on the second equation in (jl.ll) . multiplying it by 
D m d and integrating over R 3 yields 



ld_ 

Ydt 



f \D m d\ 2 dx+ f \D m+1 d\ 2 dx 
Jr 3 Jr 3 

- t D m (u-Vd)D rn d+ \D m [(\d\ 2 -l)d]D m ddx 

Jr 3 n 

[ D rn (u®d)D m+1 d+-^D m - 1 \(\d\ 2 - l)d]D m+1 ddx 



4 

Therefore 



<l I \D m+1 d\ 2 dx+ I \D m {u®d)\ 2 + C\D m - 1 [(\d\ 2 -l)d]\ 2 dx. 



(6.84) [ \D m d\ 2 dx + - f \D m+1 d\ 2 dx 

dt Jr 3 2 J R 3 

<C f \D m (u(g>d)\ 2 + \D m - 1 [(\d\ 2 -l)d]\ 2 dx. 
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The first integral on the right hand side of (|6.84|) is estimated as 

(6.85) / \D m {u® d)\ 2 dx 

<\\u\\l~\\D m d\\l 2 + \\d-w4l^\\D m u\^ 

m — 1 

+ Y,\\D i d\\l~\\D m - i u\\ 2 L2 

i=l 

<Hli-||I> m d||i a + ||d-t«o||i.||£> ni «||i a 

m—1 

+ IW^l^d-wolfft^^D™-^ 

i=l 

<\\u\\l x \\D m d\\l 2 + \\d-w \\i x \\D m u\\l 2 

m — 1 

+ j-\\D m ^d\\l 2 +C m Y / \\d-wo\\h\\D 



m — 1 

im _t oi ||2/(l-a0 
L 2 



The second integral on the right hand side of (|6.84|) is estimated as 
(6.86) 

l- Dm_1 [([d| 2 - l)d]| a da: 

m — 1 



n III -L 

/ V \D l [(d + w ){d - w )}D m - i - 1 d\ 2 dx 



(m-l)/2 



< lld-^olllooll^-^ + Cl+lld-^ollloo) y police IP^-MHi, 

4=1 

(m-l)/2 i-1 
i=l 3=1 



Using the interpolation inequality (|6.76p gives that 



(m-l)/2 

£ Wdf L A\D nM df L2 



i=l 

(m-l)/2 



<C Yl \\D m+1 d\\ 2 ^\\d-w a f L { ^ ] \\D m -^d\\ 2 L2 

(m-l)/2 



<l\\D m+i d\\i 2 + c m iM-^oiii.ii^-^ii^ 1 - 

1=1 
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Using the interpolation inequality (|6.76j) again to the last term in (|6.86j) gives that 

(m-l)/2 i_l 

E E ii^'^ii- \\ Di - jd \\l°° WD^-'dwi, 

i=l 3 = 1 

(m-l)/2 i_l 

i=i i=i 

(m-l)/2 2(2~q i ) 



1 V . 

< -p m+1 d||! 2 + c TO E N-^o||^ WD™-*- 1 * 

i— 1 

The last two inequalities combined with (|6 .86[) imply that 
(6.87) / \D m - 1 [{\d\ 2 - l)d]\ 2 dx 



^ip^dllla + lld-^olll-IP™- 1 ^!, 
(ro-l)/2 

+ (l+||d-«, ||a 8O ) ||d-t« ||i a |p m - i - 1 d||^ <1 - Ol) > 

i=l 

where we used the facts that y^f 1 > 1 an d that ||d — u>o||l 2 decays for large time. 

Finally, adding the two inequalities (|6~5T|) and ([6T51]) . and using ([6"Tg2]) . 
(j6T85|) and (|6~87|) gives that 



d 
dt 



[ \D rn u\ 2 + \D m d\ 2 dx+ [ \D m+1 u\ 2 + \D m+1 d\ 2 dx 
Jr 3 Jr 3 

< C(\\u\\ 2 Loo + \\d- wollioc + \\D 2 d\\ 2 L ^) [ \D m u\ 2 + \D m d\ 2 dx 

Jr 3 



+ C m \\d-wo\\U\\D m - 1 d\\ 2 L2 +R m , 

where 

'0, to = 1,2; 



J EZi \\D m -' l u\\ + \\d - tfldlli, \\D™~^d\\ ^ 

-H-m } 2 m-1 2 

- E™ i 1 \\d - vo\\h P m -*«ll + EST IM - ^olli, \\D™-^d\\ , 

= i?ml + Rral + Rm3 + R m 4, for 777 > 3 

with a, = 'ra^ 2 ■ The proof of the Lemma is complete. 

□ 



We are now ready to prove Theorem 12.41 by induction. 

Proof of Theorem 12.41 When 777 = 0, inequality (|2.18|) has been established 
in Lemma [4.11 Section l4l and inequality (|2.19j) is proved in Theorem 12.21 fsee [5]). 
Assume by induction that, for k = 0, 1, 2, m — 1, 

(6.88) ll^ll!=( K 3) <C K (l + t)- (k+3/2) , 

(6.89) \\D k d\\ 2 L2{R3) <C k (l+t)-^ 2 \ 
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We then apply the Fourier splitting method to the energy estimate (|6.79|) . Lemma 
4.21 Lemma \b. II and inequality (|6.78[) in Lemma EH] imply 

||«Hi- + \\d - ttfollio. + \\D 2 d\\l~ < C(l + 
Lemma 15.11 also implies 

\\d-w \\ L °a <C(l+t)-\ 
The induction hypothesis (|6.89[) implies 

WD^dWii < c m (i + t)- (m+1/2) . 

By ([57751) we have 

(6.90) ^ / \D m u\ 2 + \D m d\ 2 dx + f \D m+1 u\ 2 + \D m+1 d\ 2 dx 



^Cil + ty 1 [ \D m u\ 2 + \D m d\ 2 dx 
+ C(l + t)-^ +m +V +CR m . 



For m = 1,2, R m = 0. For m > 3, we proceed to estimate each term R m i on the 
right hand side of (|6.80p as follows. The induction hypothesis (|6.88[) and the decay 
estimate (|4.53[) in Lemma |4~T1 yield 

m/2 

(6.91) R ml =J2\\uf L2 \\D m - i u\\% {1 - ai) KCrnil+ty^-^^ 2 

< C m (l + tyi- {m+1) ^^ < C m (l + t)-(§+™+ 1 ). 

The hypothesis (|6.88[) and the decay estimate (|2.12[) in Theorem 12.21 yield 

m— l „ 

(6.92) R m3 = J2 Wd-woWhWD" 1 -^] 1 -* < C m (l + t)-- ^r— 

< c m (i + < c m (i + 

The hypothesis (|6.89[) and the decay estimate (|2.12[) in Theorem 12.21 yield 

(m-2)/2 , 



(6.93) R m2 = Y, ¥-Ml4D m - i+1 d\\v ai+1 

i=2 

3 m-i + 3/2 + 1 

<C m (l + t) 5 

< C m (l + t)-^- [m+iy ^m < C m (l + 

Similarly as to obtain (|6.93[) . we have 

(m-l)/2 _ JL _ 

(6.94) i? m4 = £ Wd-woWhWD^^dW 1 ^ <C rn {l+t)-^ +m+1 \ 

i=l 

Combing the inequalities (|6.90[) - (|6.94p yields 

(6.95) A / \D m u\ 2 + \D m d\ 2 dx + [ \D m+1 u\ 2 + \D m+1 d\ 2 dx 



dt 



< C m (l + t)- 1 f \D m u\ 2 + \D m d\ 2 dx + C m {\ + t) 
Jn 3 



-(m+5/2) 
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Applying the Fourier splitting method to the inequality (|6.95[) yields 

(6.96) / \D m u\ 2 + \D m d\ 2 dx < C m (l + t )-( m+ V 2 \ 

Jr 3 

for all m > 1. This proves inequalities (|2.18[) and (|2 . 1 9[) for all m > 1. 
Apply the following Gagliardo-Nircnberg interpolation inequality 

\\D m u\\ L ~ m < C\\D m+2 u\\l 2(M3) \\u\\l^ R3) 

witn a = 2 (m+2) • Combining (|4~55j) (in Lemma B~Tj) . (jB~W| and the last inequality 
yields 

\\D m u\\ L oe {K3) <C m (l + i)- (3/2+m+1) t-T =C m {l + t)~^. 

This proves inequality f|2 .20[) . Inequality (|2 . 2 1 1) can be proved similarly. This 
completes the proof of Theorem 12.41 
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